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Abstract 

The behavior of the specific heat c p , effective mass M*, and the thermal expansion coefficient 
a of a Fermi system located near the fermion condensation quantum phase transition (FCQPT) 
is considered. We observe the first type behavior if the system is close to FCQPT: the specific 
heat c p oc VT, M* oc while the thermal expansion coefficient a oc \[T. Thus, the 

Griineisen ratio T(T) = a/c p does not diverges. At the transition region, where the system 
passes over from the non-Fermi liquid to the Landau Fermi liquid, the ratio diverges as T(T) oc 
\j\fr. When the system becomes the Landau Fermi liquid, T(T,r) oc 1/r, with r being a 
distance from the quantum critical point. Provided the system has undergone FCQPT, the 
second type takes place: the specific heat behaves clS Cp OC VT, M* oc l/T, and a = a + bT with 
a, b being constants. Again, the Griineisen ratio diverges as T(T) oc l/x/T. 
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The nature of the strongly correlated state and the non-Fermi liquid (NFL) behavior observed in 
the high-temperature superconductors, heavy-fermion (HF) metals, two-dimensional (2D) 3 He and 
electron systems etc., are still hotly debated. It is widely believed that the observed behavior is 
determined by quantum phase transitions which occur at quantum critical points. The proximity of 
a system to quantum critical points creates its NFL behavior. The system can be driven to quantum 
critical points (QCPs) by tuning some control parameters other than temperature, such as pressure, 
magnetic field, or doping. The NFL behavior observed in many systems suggests that in principle 
there can exist many CQPs. The situation becomes more misty and obscure, if one takes into account 
that when a system is close to QCP we are dealing with the strong coupling limit where an absolutely 
reliable answer cannot be given, based on purely first principle theoretical grounds. Therefore, the 
only way to verify what type of QCP occurs is to consider experimental facts which describe the 
behavior of a system. Only recently, experimental facts have appeared which deliver experimental 
grounds to understand the nature of quantum phase transition producing the NFL behavior of HF 
metals. 

Recent measurements for the HF metal CeNi2Ge2 have shown that the thermal expansion coeffi- 
cient Oi(T) exhibits the NFL behavior, oc(T) = a^/T+bT, the specific heat c p (T)/T = 7 — cVT+d/T 3 , 
and that the Griineisen ratio T(T) = a(T)/c p (T) diverges as T(T) oc 1/y/T PQ. Here, a, 6, 70, c, d are 
constants, and the low-temperature upturn in the specific heat related to the d/T 3 term can be as- 
cribed to the high-temperature tail of a Shotky anomaly, whose influence on the measured Griineisen 
ratio is small pQ. While measurements on the HF metal YbRh2(Sio.9sGeo.o5)2 have revealed the 
different NFL behavior: at T < 1 K, the thermal expansion coefficient a(T) = a\ + aoT, and the 
Griineisen ratio diverges as T(T) oc l/T 2 , where the exponent z is fractional, z = 0.7 ± 0.1 pQ. 

A general consideration based on the assumption of scaling and a model study on the spin-density- 
wave j2] has produced an explanation for the NFL behavior of the thermal expansion coefficient in 
the case of CeNi 2 Ge2 but failed to explain the NFL behavior observed in YbRh 2 (Sio.95Geo.o5)2- As a 
result, a challenge still remains how to explain within one conceptual framework the NFL behavior 
observed in these HF metals, both of which are considered to be close to QCPs. 

In this Letter, we show that within the framework of the fermion condensation quantum phase 
transition (FCQPT) [3IIj, which leads to the formation of fermion condensate (FC) [5], it is possible 
to understand the NFL behavior observed in different strongly and highly correlated Fermi systems 
such as high-T c superconductors [31111310], heavy-fermion metals jHEllH], and 2D systems [USHII]. 
We consider the specific heat c p , effective mass M*, and the thermal expansion coefficient a of 
a Fermi system near FCQPT and show that at sufficiently high temperatures, the specific heat 
c p oc VT, M* oc 1/yT, and the thermal expansion coefficient a(T) oc \/T. Thus, the Griineisen 
ratio T(T) = a/c p does not diverge. At the transition temperatures, when the system comes from 
the NFL behavior to the Landau Fermi liquid (LFL) one, the ratio diverges as T(T) oc 1/VT. As 
the system possesses the LFL behavior at low temperatures, T(T,r) oc 1/r, here r is a distance 
from FCQPT. If the system has undergone FCQPT and FC is formed, the specific heat behaves 
as Cp oc y/T, and a = a + bT with a, b are constants. Again, the Griineisen ratio diverges as 
T(T) oc 1/Vr. Thus, the second type of behavior can be observed in a Fermi system with FC, which 
nonetheless demonstrates the main features of the critical behavior, as if the system were located 
near QCP. Our results are in good agreement with measurements on single crystals of CeNi 2 Ge2 and 
YbRh 2 (Sio.95Ge .o5)2 PP- We predict that if the system is driven back to LFL by the application 
of a magnetic field B, then in the case of YbRh 2 (Sio.95Geo.o5)2 the thermal expansion coefficient 
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a(B) oc 1/y/B, while in the case of CeNi 2 Ge2 the coefficient behaves as a(B) oc B~ 2 ^ 3 . 

Let us start by considering the key points of the FC theory. We concentrate on the 3D case 
because heavy fermion metals are 3D structures ^U]. On the other hand, there is no a qualitative 
difference between 3D and 2D cases in the theory of FC, and our results are applicable to both 2D 
and 3D cases [3 |H] . FC represents a new solution of the Fermi liquid theory equations [TT] for the 
quasiparticle distribution function n(p, T) 

S(F-fiN) l-rc(p,T) 

= e p,T )-fj, T -Tin — - — = 0, (1) 

<m(p,T) n(p,T) 

which depends on the momentum p and temperature T. Here F = E — TS is the free energy, and 
\i is the chemical potential, while e(p, T), 



£ (p,T) = |f% (2) 
5n(p,T) 

is the quasiparticle energy. This energy is a functional of n(p, T) just like the total energy E[n(p)], 
entropy S^nQo)] and the other thermodynamic functions. The entropy 5*[n(p)] is given by the familiar 
expression 

S[n(p)} = -2 J [n(p, T) In n(p, T) + (1 - n(p, T)) ln(l - n(p, T))] (3) 

which stems from purely combinatorial considerations. Eq. (1) is usually presented as the Fermi- 
Dirac distribution 

Xp,t)-//) 



n(p,T) = il +exp 



T 



(4) 



At T — ► 0, one gets from Eqs. (1) and (3) the standard solution np(p, T — ► 0) — ► ^(^f — p), with 
e(p ~ pi?) — /i = Pf(p — Pf)/MI, where is the Fermi momentum, #(pf — V) is f ne step function, 
and Ml is the Landau effective mass [TT] 

1 _ l &(p,T->0) 

It is implied that in the case of LFL M£ is positive and finite at the Fermi momentum pp. As a 
result, the T-dependent corrections to M£, to the quasiparticle energy e(p), and other quantities, 
start with T 2 -terms. But this solution is not the only one possible. There exist solutions of Eq. (1) 
associated with the so-called fermion condensation 0. Being continuous and satisfying the inequality 
< rio(p) < 1 within some region in p, such solutions ^o(p) admit a finite limit for the logarithm in 
Eq. (1) at T — > yielding 

e(p)-(i = 0, if 0<n (p) <l;Pi<P<Pf ■ (6) 

At T — > 0, Eq. (6) defines a new state of electron liquid with FC lEj which is characterized by 
a flat part of the spectrum in the (pf — Pi) region and can strongly influence measurable quantities 
up to temperatures T <C Tf [3 E]- Note, that a formation of the flat part of the spectrum has 
been recently confirmed ^3 E], and the momenta pi and pf have to satisfy pi < pp < Pf- At 
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T — > 0, Eq. (6) defines a particular state of a Fermi liquid with FC, for which the order parameter 



k (p) = y (1 — n o(p))^o(p) has finite values in the (p/ — Pi) region, whereas the superconducting gap 
Ai — > in this region. Such a state can be considered as superconducting, with an infinitely small 
value of Ai, so that the entropy S(T = 0) of this state is equal to zero. It is obvious that this state 
being driven by the quantum phase transition disappears at T > [3] . When the density x tends to 
QCP located at Xpc x ~^ x fc, the Landau amplitude F L (p = Pf,Pi = Pf,%) as a function of the 
density x becomes smaller, and the flat part vanishes. Then, at T — > 0, Eq. (6) has the only trivial 
solution e(p = pf) = and the quasiparticle occupation numbers are given by the step function, 
n(p) = 6{p F -p) 0. 

It follows from Eq. (6) that at Tf ^> T > 0, the system becomes divided into two quasiparticle 
subsystems: the first subsystem is occupied by normal quasiparticles with the finite effective mass M F 
independent of T at momenta p < p iy and the second subsystem in the (pf—pi) range is characterized 
by the quasiparticles with the effective mass M FC (T) [3| IT5] 

M FC c PF P±^. (7) 

There is the energy scale Eq separating the slow dispersing low energy part, related to the effective 
mass M FC , from the faster dispersing relatively high energy part, defined by the effective mass M£. 
It follows from Eq. (7) that Eq is of the form [B] 

Eq ~ AT. (8) 

The described system can be viewed as the strongly correlated system and has the second type of the 
behavior (the first type will be considered below). At T — > 0, the system demonstrates the typical 
critical behavior as if the system were close to a quantum critical point: the effective mass diverges 
as M* oc 1/T, while the application of magnetic field B restores the LFL behavior with the effective 
mass M* oc \j\l~B 0. Note that such a behavior is in agreement with experimental facts [16J. 

The existence of the two quasiparticle subsystems can be illuminated by calculating the thermal 
expansion coefficient u(T) [T^], which is given by ^7] 




x ( d(S/x 
3K { dx 



I • (9) 

T 



Here, P is the pressure and V is the volume. The compressibility K is not expected to be singular 
at FCQPT and in systems with FC, because FC is attached to the Fermi level, and it moves along 
as [a(x) changes, while the compressibility K = dp,/d(Vx) is approximately constant JH]- Inserting 
Eq. (3) into Eq. (9), we find that 

S(T) 1 f dnfaT) ^ l-n(p,T) dp 
a(T) =3Kx--3kJ -^x~ ln n(p,T) (2^ 

Because the region (pf—pi) is occupied by FC, we have from Eq. (1) that the ratio (3 = (e(p)— fj)/T ~ 
ln(l — no(p))/ri(p) does not depend on the temperature. We expect corrections to be of the order 
T at Tf ^> T, so that (3 ~ const + T. Thus, the right hand side of Eq. (10) can be estimated as 
oc a + bT, the term a is determined by the FC contribution, and bT is given by normal quasiparticles 
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with the effective mass M£. Therefore, in the case of systems with FC, the thermal expansion 
coefficient otpc — a + bT, while in the case of LFL, one obtains from Eq. (10) the standard expression 
a{T) ~ MIT/ppK. On the other hand, for the two quasiparticle subsystems with M FC oc 1/T, see 
Eq. (7), Eq. (10) translates into 

We again find that apc{T) = a + bT. This behavior of the thermal expansion coefficient predicted 
in Ref. [15J is in good agreement with experimental measurements on YbRh 2 (Si .95Ge .o5)2 PQ- 

The application of magnetic fields to YbRh 2 (Sio.95Ge .o5)2 restores the LFL behavior, and the 
effective mass scales as M* cx 1/yB at sufficiently low temperatures T < T*(B), where T*(B) oc \J~B 
is the temperature above which the system comes back to the NFL behavior [7j. As a result, we have 
from Eq. (10) that the thermal expansion coefficient depends on the magnetic field as ot(B) oc 1/yB. 

According to the Landau theory, the specific heat per unit volume is 

c(rH2/(e(p) _^_|L (12) 

At low temperatures, as in the common case, we can ignore the difference 5c between c v and c p in 
systems with FC because (c„(T) — c p {T)) oc Ta 2 (T) [T7], and as it follows from the above, 5c oc T [15j. 
On the other hand, this difference is large when compared to the common case for 5c oc T 3 |17j . As 
we have seen above, the ratio (3 is approximately constant in the region (pf —pi). Using this fact and 
differentiating both sides of Eq. (4) with respect to T, we can check that dn(p)/dT oc T. As a result, 
it is seen from Eq. (12) that the FC region does not contribute to the singular behavior of the specific 
heat. The contribution coming from the normal quasiparticles with M£ is of the standard form being 
proportional to T. There is a specific contribution related to the spectrum e(p) which insures the 
connection between the dispersionless region (pf —p^ occupied by FC and the normal quasiparticles 
located at p < p^ and at p > pf, and therefore it is of the form e(p) ~ (p — pf) 2 ~ {pi — p) 2 . 
Such a form of the spectrum can be verified in exactly solvable models for systems with FC [5]. 
One can check, that the contribution of this spectrum to c p oc a/T, and this consequence coincides 
with one obtained in exactly solvable models 0. As a result, we have c p (T) oc ai^/T + b\T, with 
a\,b\ being some constants. It is worth noting that when calculating quantities like c p containing 
the derivative dn(p)/dT, it is incorrect to insert Eq. (7) into the standard formula for the specific 
heat, c(T) oc M*T. Now we can calculate that the Griineisen ratio diverges as T(T) oc 1/yT. 
The obtained Griineisen exponent z = 1/2 is close to the experimental one z = 0.7 ± 0.1 1. We 
conclude that in contrast to the common belief, the heavy-fermion metal YbRh 2 (Sio.95Geo.o5)2 need 
not be located very near CQP, and there is no need to introduce a new type of QCP to explain its 
properties. 

Consider the case when the system is close to FQCPT and demonstrates the first type of the 
behavior. FCQPT can be induced by tuning the density x. When the system approaches FCQPT 
from disordered phase, FCQPT manifests itself in the divergence of the quasiparticle effective mass 
M* as the density tends to the critical density xpc, or the distance r ~ \x — xpc\ —> [El HI] 



Note, that Eq. (13) is valid in both cases of 2D and 3D systems |S]. Since the effective mass M* is 
finite, the system exhibits the LFL behavior at low temperatures T ~ T*(x) oc \x — xpc\ 2 |S|- At 
sufficiently high temperatures, the system possesses the NFL behavior. At \x — xpc\l x fc <S 1> this 
behavior can be viewed as the highly correlated one, because the effective mass strongly depends on 
the control parameters such as the density, temperature and magnetic fields .$]. In contrast to the 
strongly correlated system, the highly correlated system does not have the energy scale given by Eq. 

(8) , Eq. (7) is not also valid, and at T < T*(x) it behaves as LFL [Hj. Let us calculate the Griineisen 
ratio at T ~ T*(x). The entropy S is given by the LFL formula, S oc M*T. Inserting this into Eq. 

(9) and taking into account Eq. (13), we obtain a(T) oc T/\x — xpc\ 2 while c p oc T/\x — Xpc\- As a 
result, the Griineisen ratio diverges as 

T(T, r) oc r oc - oc M*. (14) 

\x — Xfc\ r 

At the transition temperatures, the system passes from the NFL behavior to the LFL one, and the 
effective mass depends on both T and x. In this case, we expect Eq. (14) to hold true because the 
effective mass has to be continuous over the transition region. 

Now consider the Griineisen ratio at T ^> T*(x) when the effective mass depends predominantly 
on the temperature. Landau equation relating the mass M of an electron to the effective mass of 
the quasiparticles is of the form ^T] 

W* = M + I Fl(p ' Pl ' ^ (^3- ( 15 ) 
Applying Eq. (15) at T < T*(x), we obtain the common result 

M " = l-Ji M /3 - < 16 ' 

Here N is the density of states of a free Fermi gas and Fl(x) is the p-wave component of the Landau 
interaction. At x — > x F c, the denominator in Eq. (16) tends to zero and one obtains Eq. (13). The 
temperature smoothing of the step function 0(p F —p) at p F = (x/37r 2 ) 1//3 induces the variation of the 
Fermi momentum Ap F ~ TM*/pp. We assume that the amplitude Fl has a short range go ^ Pf of 
interaction in the momentum space, which is a common condition leading to the existence of FC and 
nearly-localized Fermi liquids jHHSOl- If the radius is that q ~ Ap F ~ T M*/p F , corrections to the 
effective mass are proportional to T at T ~ T . Here T oc \x — x F c\ is a characteristic temperature 
at which the system's behavior is the NFL one. On the other hand, since T*(x) To, we have 
go ^> T*(x)M* /pf, and the system behaves as LFL at T ~ T*(x), that is corrections to the effective 
mass start with T 2 terms. We can also conclude that the transition region is rather large compared 
with T*(x), being proportional to T . 

In the case of T ~ T , we again can use Eq. (16) with Fl(p F + Ap F ) ~ Fl(p F ) + AAp F , where 
A oc dF F (x)/dx. Substituting this expansion of Fl(p F + Ap F ) into Eq. (16) we find that 

MM 

AAp F M*T v ' 
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In deriving Eq. (17) it is implied that the system is close to FCQP so that (1— N Fl(p F )) <C N AAp F . 
We can say that at T ~ T , Ap F induced by T becomes larger then the distance r from FCQP, and 
Ap F > \pp C — Pf\, where p F c corresponds to x F c- As a result, Eq. (17) becomes 

M* oc -=. (18) 
VT 

Then from Eqs. (10) and (18), we obtain that at T ~ T , the main contribution to the thermal 
expansion coefficient a(T) oc a\VT + b±T in accordance with experimental facts pQ. Here the b\ term 
comes from the contribution of normal quasiparticles. 

By applying magnetic fields the system in consideration can be driven back to LFL with the 
effective mass M* oc B~ 2 / 3 at temperatures T < T*(B), where T*(B) oc 5 4/3 jH]. Thus, in the 
case of CeNi 2 Ge2, it follows from Eq. (10) that the thermal expansion coefficient depends on the 
magnetic field as a(B) oc B~ 2 / 3 . We note that the magnetic field should not be too high, otherwise 
the distance r becomes large, and the system cannot be considered to be close to the quantum critical 
point. A more detailed analysis of this issue will be published elsewhere. 

When considering the specific heat, we again ignore the difference (c p (T) — c v (T)) oc Ta 2 (T) oc T 2 . 
Note that this difference is large in comparison to the common case. It follows from Eq. (12) that 
the specific heat c p (T) ~ c v (T) oc y/T. This result implies that the Griineisen ratio does not diverge 
at T ~ To. However, at the transition region where the effective mass depends on both T and x, 
it is seen from Eq. (14) that the ratio diverges, T(T) oc 1/y/T, and finally at T ~ T*(x), we have 
T(T, r) oc l/\x — Xpc I ^ l/ r - 

It follows from measurements on CeNi 2 Ge 2 that T(T) oc 1/VT, while the behavior of the specific 
heat c p (T) oc 70T — cT 3 / 2 pQ suggests that the system in question enters the transition region. Both 
these observations are consistent with our consideration. Therefore, we expect that at decreasing 
temperatures the ratio will be saturated. 

In conclusion, we have shown that our simple model based on FCQPT explains the critical 
behavior of the thermal expansion coefficient, specific heat and the Griineisen ratio observed in 
different heavy fermion metals. In the case of CeNi 2 Ge 2 , the behavior of these quantities can be 
explained by its closeness to FCQPT, while the electronic system behaves as a highly correlated 
liquid. In the case of YbRh 2 (Sio.95Geo.o5)2 the behavior of these quantities is different, and can be 
explained by the presence of FC in this metal, that is the electronic system has undergone FCQPT 
and behaves as a strongly correlated liquid. We also predict that if the system is driven back to LFL 
by the application of a magnetic field B, then in the case of a strongly correlated liquid the thermal 
expansion coefficient a(B) oc l/v^B, while in the case of a highly correlated liquid the coefficient 
behaves as a(B) oc B~ 2 / 3 . 
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